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ABSTRACT 

We investigate what may be the origin of the presently observed spatial distribution of 
the mass of the Galactic Old Halo globular cluster system. We propose its radial mass 
density profile to be a relic of the distribution of the cold baryonic material in the 
protoGalaxy. Assuming that this one arises from the profile of the whole protoGalaxy 
minus the contribution of the dark matter (and a small contribution of the hot gas 
by which the protoglobular clouds were bound), we show that the mass distributions 
around the Galactic centre of this cold gas and of the Old Halo agree satisfactorily. In 
order to demonstrate our hypothesis even more conclusively, we simulate the evolution 
with time, up to an age of 15Gyr, of a putative globular cluster system whose initial 
mass distribution in the Galactic halo follows the profile of the cold protogalactic gas. 
We show that beyond a galactocentric distance of order 2 to 3 kpc, the initial shape 
of such a mass density profile is preserved inspite of the complete destruction of some 
globular clusters and the partial evaporation of some others. This result is almost 
independent of the choice of the initial mass function for the globular clusters, which 
is still ill-determined. The shape of these evolved cluster system mass density profiles 
also agree with the presently observed profile of the Old Halo globular cluster system, 
thus strengthening our hypothesis. Our result might suggest that the flattening shown 
by the Old Halo mass density profile at short distance from the Galactic centre is, at 
least partly, of primordial origin. 
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1 INTRODUCTION 

Globular Clusters (GC) are thought to be the oldest bound 
stellar systems in our Galaxy. Their study provides there- 
fore valuable information about the early Galactic evolution. 
In this respect, a major problem is that we do not know 
whether what we presently observe is still representative of 
the initial conditions and, thus, a fossil imprint of the for- 
mation process, or whether the initial conditions have been 
wiped out by a 15 Gyr long evolution within the tidal fields 
of the Milky Way. Modelling the dynamical evolution of the 
Galactic Globular Cluster System (CCS) is thus of great in- 
terest as it helps us to go back in time to the earliest stages 
of the cluster system and to disentangle the formation and 
evolutionary fingerprints (see, e.g., Gnedin & Ostriker 1997, 
Baumgardt 1998, Vesperini 1998, Fall & Zhang 2001). The 
GCs most vulnerable to evaporation and disruption are the 
low-mass clusters located at small galactocentric distance. 
As a result, the evolution with time of a GCS is markedly 
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determined by the initial distribution of the GCs in space 
around the Galactic centre as well as by their initial mass 
spectrum. 

As for the presently observed spatial distribution of the 
Galactic halo GCs, it is centrally concentrated with the den- 
sity varying as D~ 3,5 (D is the Galactocentric distance) over 
most of the halo (Zinn 1985). In the inner 3-4 kpc, the distri- 
bution flattens to something closer to an D~ 2 dependence. 
As a result, the overall distribution is conveniently described 
by a power-law with a core (see Section 2). The observed cen- 
tral flattening probably arises from a combination of several 
effects: our failure to discover some GCs in the heavily ab- 
sorbed central regions of the Galaxy, distance errors, and the 
real flattening of the distribution. It is still unclear whether 
such a flattening is of primordial origin and reflects the ini- 
tial spatial distribution of the system, or whether it has been 
completely determined by evolutionary processes. The lat- 
ter are especially effective at small galactocentric distances 
where the GC relaxation time is short, causing the disrup- 
tion of some GCs and the partial evaporation of some others. 
More generally, as far as we are aware of, the shape of the 
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whole spatial distribution of the halo GCS has received no 
explanation. 

The aim of this paper is to present a possible model ac- 
counting for the distribution around the Galactic centre of 
the mass content of the halo GCS, that is, the variation with 
galactocentric distance of the mass density of the halo GCS. 
We will consider solely the Old Halo (OH) GCs, excluding 
from our analysis the Younger Halo (YH) subsystem. This 
one is made of GCs suspected of having been accreted and 
is thus of limited relevance to the earliest stages of the main 
body of our Galaxy (see, e.g., Zinn 1993, van den Bergh 
1993). A common trend of different scenarios for the for- 
mation of GCs is to assume that their gaseous progenitor 
clouds, whatever their size (~ 1O 6 M , e.g., Fall & Rees 
1985; ~ 10 9 M Q , e.g., Harris & Pudritz 1994), are embed- 
ded in a hot and tenuous background at the virial temper- 
ature. We build on this picture, splitting the protoGalaxy 
in a set of three components: a gaseous medium consisting 
of a hot and a cold phases (in the form of a collection of 
cold and dense clouds, the gaseous precursors of the halo 
GCs, pressure-bound by a hot medium) , and a dark matter 
corona. We investigate whether the mass density profile of 
the halo GCS may trace the profile of the cold phase, that 
is, of the cold baryonic matter which was available to star 
formation some 15 Gyr ago. 

The outline of the paper is as follows. In Section 2, 
we build the radial mass density profile of the OH and ob- 
tain new fits of a power-law with a core, with parameter 
values appropriate to this subsystem of GCs only. We also 
summarize the evidence from the literature following which 
mass-related quantities may be better probes to the initial 
conditions than number related quantities. In Section 3, we 
present our hypothesis regarding the shape of the initial 
mass density profile of the OH GCS and we compare the OH 
profile obtained in Section 2 with our suggested model, that 
is, with the radial density profile of the cold protogalactic 
material. In Section 4, we simulate the evolution with time 
of a GCS whose initial mass density profile mirrors the one 
of the cold protogalactic gas, and this for various initial GC 
mass spectra. We compare the spatial mass distributions of 
these evolved systems with the presently observed spatial 
mass distribution of the OH. Finally, our conclusions are 
presented in Section 5. 



2 RADIAL DISTRIBUTION OF THE HALO 

GCS MASS 

The observed radial distribution of the Galactic GCS (i.e., 
the number of GCs per unit volume in space as a function 
of Galactocentric distance D) is often parametrized by a 
simple power law with a core: 

P ( D) = P0 (i + §-y\ a) 

Realizing a "chi-by-eye" fitting of this type of curve on the 
observed spatial distribution of the GCS, Djorgovski & Mey- 
lan (1994) found good matches for 7 ~ 3.5-4.0 and D c ~ 0.5- 
2kpc, the steepest slope being associated with the largest 
core. This approach is purely empirical however and is not 



meant to imply any physical meaning of the distribution 
given by equation (1). 

Both the number and the mass density profiles of the 
halo ([Fe/H] < -0.8) GCS can be described by equation (1), 
as both profiles are indistinguishable in shape. However, as 
argued by McLaughlin (1999), the spatial distribution of the 
mass is likely to be a better estimate of the initial conditions 
than its number counterpart. Dynamical evolution target- 
ting mostly low-mass clusters and these ones accounting for 
a limited fraction of the GCS mass (see below and Section 
4), the decrease with time of the total GCS mass is much 
slower than the decrease of the total number of GCs. 

This interesting property comes from how the shape of 
the initial mass spectrum of the halo GCs may have looked 
like. In our Galaxy, the luminosity function of the halo GCs 
1 (the number of GCs per unit absolute magnitude, which is 
proportional to the number of objects per logarithmic mass 
interval) is bell-shaped and usually fitted with a gaussian. 
However, the underlying mass spectrum (i.e., the number 
of objects per linear mass interval) is well fitted by a two- 
index power-law, with exponents ~ —2 and ~ —0.2 above 
and below ~ 1.5 x 10 5 M , respectively (McLaughlin 1994). 
The peak of the gaussian magnitude function in fact coin- 
cides with the cluster mass at which the slope of the mass 
spectrum changes. The slope of the high mass regime is 
reminiscent of what is observed in interacting and merg- 
ing galaxies (see, e.g., Whitmore & Schweizer 1995, Whit- 
more et al. 2002) where systems of young GCs show well 
defined power-law with slopes ranging between -1.8 and -2 
for their luminosity spectrum (but see the discussion in Sec- 
tion 4). This thus suggests that the initial mass spectrum 
of the halo GCS may have been itself a single-power law. 
Numerous studies of GCS dynamical evolution modelling 
have shown that a Hubble time long evolution turns such 
an initial spectrum into the presently observed one (e.g., 
Baumgardt 1998, Vesperini 1998, Fall & Zhang 2001). In 
fact, low-mass clusters being the most vulnerable to evapo- 
ration and disruption, the GC mass spectrum gets severly 
depleted below a turnover of ~ 1.5 x 10 5 M , leading to a 
much shallower mass spectrum (i.e., slope ~ —0.2) in the 
low-mass regime. Pal 5 constitutes a striking example of a 
low-mass GC currently dissolved by the Galactic tidal fields 
(Odenkirchen et al., 2001; Dehnen et al., 2004). 

Assuming that most of the GCs more massive than the 
turnover are spared by the dynamical evolution, McLaugh- 
lin (1999) compares the GCS initial and final mass spectra 
(i.e., the single power-law with the two- index power-law) and 
derives useful formulae to estimate the fraction of surviving 
clusters, both in term of mass and numbers (his equations 4- 
7) . His results show that mass-related quantities are reason- 
ably preserved by a Hubble time long evolution, even though 
low-mass GCs are disrupted in large numbers. Considering 
the specific case of the Milky Way (power-law slopes of -0.2 
and of around -1.8 to -2.0 for the low and high-mass regimes 
respectively, a turn-off mass of 1.5 x 1O 5 M , and lower and 

1 In what follows, we adopt the nomenclature of McLaughlin 
& Pudritz (1996). We call mass/luminosity spectrum the num- 
ber of objects per linear luminosity /mass interval, dN/dm or 
dN/dL, while we refer to the mass/luminosity function to de- 
scribe the number of objects per logarithmic luminosity/mass in- 
terval, dN/d log m or dN/d log L. 
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Table 1. OH sample 



NGC288 


NGC5286 


NGC5986 


NGC6218 


NGC6293 


HP1 


NGC6522 


NGC6638 


NGC6752 


Pal2 


NGC5466 


NGC6093 


NGC6235 


NGC6341 


NGC6362 


NGC6535 


NGC6652 


NGC6779 


NGC1904 


NGC5634 


NGC6121 


NGC6254 


NGC6325 


NGC6402 


NGC6540 


NGC6656 


NGC6809 


NGC2298 


NGC5694 


NGC6101 


Pall5 


NGC6333 


NGC6401 


NGC6544 


NGC6681 


NGC6864 


NGC4372 


NGC5824 


NGC6144 


NGC6266 


NGC6355 


NGC6397 


NGC6541 


NGC6712 


NGC7078 


NGC4833 


NGC5897 


NGC6139 


NGC6273 


IC1257 


NGC6426 


NGC6558 


NGC6717 


NGC7089 


NGC5024 


NGC5904 


NGC6171 


NGC6284 


NGC6366 


NGC6453 


NGC6569 


NGC6723 


NGC7099 


NGC5053 


NGC5946 


NGC6205 


NGC6287 


Ter4 


NGC6517 


NGC6626 


NGC6749 


NGC7492 



upper mass limits of 1O 4 M and 10 6 Mq, respectively), the 
application of his formulae shows that the initial mass of the 
GCS is decreased by 40 per cent, while the fraction of sur- 
viving clusters is 16 per cent. In case of a lower limit for the 
cluster initial mass range, say 1O 3 M0, the contrast between 
the decreases in mass and number is even more striking, i.e., 
the survivors still represent 44 per cent of the initial total 
mass but 2 per cent of the initial number only. McLaughlin 
(1999) 's formulae thus illustrate that, compared to number 
related quantities, mass related quantities are less markedly 
affected by a 15 Gyr long evolution. As a result, the radial 
mass density profile of the halo GCS may be considered as 
a reasonably reliable estimator of their initial distribution 
around the Galactic centre. We will make this point more 
quantitative in Section 4 and show that such an hypothesis 
is indeed robust. 

In this paper, we are interested in understanding the 
origin of the initial spatial distribution of the halo GCS 
mass within the Milky Way, which we approximate by the 
presently observed mass density profile at this stage of the 
discussion. We do not consider the more metal-rich, presum- 
ably second generation, bulge/disc GCs ([Fe/H] > —0.8). 
Also, the halo ([Fe/H] < —0.8) subsystem itself could be 
divided into two groups, traditionally referred to as the OH 
and the YH (see, e.g., Zinn 1993, Van den Bergh 1993, 
Mackey & Gilmore 2004) . Evidence supporting the existence 
of such two distinct halo subsystems have been accumulat- 
ing over the past years. OH and YH GCs show differences in 
horizontal branch morphology, age, kinematics, spatial dis- 
tribution (see Parmentier et al. 2000, their Section 2, for 
a review), as well as differences in the distribution of their 
core radius (Mackey & Gilmore 2004) . The properties of the 
OH group are consistent with the majority of its members 
having been formed "in situ" during the large-scale collapse 
of the protogalactic cloud, as envisioned by Eggen, Lynden- 
Bell & Sandage (1962). On the other hand, the YH GCs 
are not native to the Galaxy, probably having been formed 
in external dwarf galaxies and afterwards accreted into the 
outer halo while their host galaxies were being swallowed 
by the Milky Way, as suggested by Searle & Zinn (1978). 
The current accretion of the Sagittarius dwarf galaxy and 
of its small GCS is the smoking gun of this process. As we 
are interested in the mass density profile of the GCS which 
formed within the original potential well of the Galaxy, we 
restrict our attention to the OH GCs. This OH/YH divi- 
sion has already been proven most fruitful as the OH GCS 
shows a metallicity gradient and obeys a mass-metallicity 
relation, features predicted by simple self-enrichment mod- 



els (Parmentier et al. 2000 and Parmentier & Gilmore 2001, 
respectively), while the whole halo GCS (OH+YH) does not. 

Lists of OH and YH GCs are provided in Lee et 
al. (1994) and Da Costa & Armandroff (1995). With re- 
spect to these, we have made two slight changes however. In 
our present OH sample, we ignore NGC 2419. Although co- 
eval with the inner halo (Harris et al., 1997; Salaris & weiss 
2002), this GC is located at a galactocentric distance of or- 
der 90 kpc and is thus unlikely to belong to the main body 
of the Galaxy. Moreover, van den Bergh & Mackey (2004) 
show that NGC 2419 and u> Cen on the one hand, and the 
other halo GCs on the other hand, are at different loeii in a 
half-light radius vs absolute visual magnitude diagram. They 
thus suggest that, as u> Cen (which we also exclude from our 
sample), NGC 2419 might be the tidally stripped core of a 
former dwarf spheroidal galaxy. An efficient tidal stripping 
would however require NGC 2419 to cross the inner Galactic 
regions. Unfortunately, its orbit is still ill-determined. Also, 
unlike lo Cen, there is no evidence for a metallicity spread 
among the cluster giants. If NGC 2419 is actually the rem- 
nant of a former dwarf galaxy, the parent galaxy might, like 
the Ursa Minor dwarf galaxy (van den Bergh 2000), have 
produced a single generation of metal-poor stars. Bearing 
these caveats in mind, we thus note that the main pecu- 
liarity of NGC 2419 with respect to the bulk of the OH is 
its large galactocentric distance. Neglecting NGC 2419, the 
OH is thoroughly contained within D<40 kpc. Additionally, 
we have moved the cluster NGC 6864 from the YH group 
to the OH group. According to the former (1999) edition 
of the GC McMaster Catalog (Harris 1996), the metallicity 
and the horizontal branch ratio (HBR) of NGC 6864 are 
-1.32 and -0.42, respectively. The updated (2003) values 
being [Fe/H] =-1.16 and HBR=-0.07, the location of this 
cluster in the [Fe/H] vs HBR diagram shows that it is more 
likely a member of the OH group rather than a YH GC. Fi- 
nally, about 25 halo GCs of the Harris Catalogue still miss 
HB index measurements. Using recently published colour- 
magnitude diagrams, (e.g., Piotto et al., 2002), Mackey & 
Gilmore (2004) have sorted most of these yet undefined GCs. 
All the GCs to which Mackey & Gilmore (2004) have as- 
signed an OH membership have been added to our OH sam- 
ple. This one is presented in Table 1. 

Before proceeding further, we derive estimates for the 
parameters 7 and D c appropriate for the sole OH, as previ- 
ous fits refer to either the whole Galactic GCS or the whole 
halo system, thus ignoring their heterogeneity. We fit 

Log w p{D) = Logw po ~ 7 Log 10 (l + 77- ) (2) 
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Table 2. Results of fitting a power-law with a core on the ob- 
served mass density profile of the OH subsystem. 
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Figure 1. Fits of a power-law with a core on the observed mass 
density profile of the OH subsystem. The different curves corre- 
spond to the values of the parameters as given in Table 2 



to the observed OH mass density distribution through a 
Levenberg-Marquardt algorithm (Press et al. 1992). Our 
source for the galactocentric distances D and the absolute 
visual magnitudes M„ is the McMaster database compiled 
and maintained by Harris (1996, updated February 2003). 
Cluster absolute visual magnitudes have been turned into lu- 
minous mass estimates by assuming a constant mass-to- light 
ratio m/L„=2.35. This value corresponds to the average of 
the mass-to-light ratios of the halo GCs for which Pryor & 
Meylan (1993) derived dynamical mass estimates. The mean 
m/L„ is almost independent of the halo group considered, 
that is, OH or OH+YH. 

The OH mass profile is derived by binning the data with 
two different bin sizes: AlogD = 0.1 and 0.2 (D is in kpc), 
corresponding to 16 and 8 points, respectively. As for the size 
of the error bars, a Poissonian error on the number of GCs 
in each bin is combined with a fixed error on the mass-to- 
light ratio. In fact, not all GCs show the same mass-to-light 
ratios, the standard deviation in the Pryor & Meylan (1993) 
compilation being of order <x Log (M/z,„) = 0.17. 

Results of the fitting procedure are presented in Table 
2 and superposed to the observed distribution in Fig. 1. For 
each fit, we also give the \ 2 an d the incomplete gamma func- 
tion Q(i//2,x 2 /2) (y is the number of degrees of freedom) 
which provides a quantitative measure for the goodness-of- 
fit of the model (Press et al. 1992). Keeping all three pa- 
rameters free, we obtain a slope much steeper (7 — —5) than 
suggested by previous works (7 ~ —3.5 to — 4.0; Zinn 1985, 
Djorgovski & Meylan 1994). 

Our slope is coupled with a rather large core (D c ~ 
3.6 kpc), however. In fact, functions like equation (1) or (2) 
show a core-slope degeneracy in the sense that two distinct 
fits can provide a satisfactory agreement with the data pro- 
vided that a steeper slope is associated with a larger core. 
This is also reflected by the large error bars for D c and 7. 
Therefore, when describing a density profile through equa- 

2 We remind the reader that a Q value of 0.1 or larger indicates 
a satisfactory agreement between the model and the data 



ALogD 


Logio(po) 


D c 


-7 


x 2 


Q 


0.1 


5.6 ± 0.3 


3.7 ± 2.0 


-5.1 ± 0.9 


11.7 


0.55 


0.2 


5.6 ± 0.4 


3.5 ± 2.4 


-4.9 ± 1.1 


3.8 


0.58 


7 imposed 


0.1 


6.4 ± 0.4 


0.9 ± 0.3 


-3.5 


17.7 


0.22 


0.2 


6.3 ± 0.5 


1.0 ± 0.4 


-3.5 


7.0 


0.32 


0.1 


6.0 ± 0.3 


1.7 ± 0.4 


-4.0 


13.9 


0.46 


0.2 


5.9 ± 0.3 


1.9 ± 0.5 


-4.0 


4.8 


0.57 



tion (1) or (2), it is important to quote the slope as well as 
the core length. We have thus performed additional fits in 
which the value of the exponent is set to values previously 
quoted in the literature, i.e., -3.5 and -4. Such shallower 
slopes again provide satisfactory agreement with the data, 
as indicated in Table 2. With respect to Djorgovski & Mey- 
lan (1994), who fitted the whole Galactic GCS (i.e., disc plus 
halo), core sizes are not markedly different. 



3 ORIGIN OF THE OH GCS MASS DENSITY 
PROFILE 

As already mentioned, the description of the OH GCS by 
equation (1) or (2) is purely empirical and does not imply 
any physical meaning. As far as we are aware of, no explana- 
tion has been put forward to explain the shape of the radial 
mass density profile poh{D). In this paper, we propose that 
Poh(D) is tracing the mass distribution of the gas which, 
some 15 Gyr ago, was available to the process of star forma- 
tion. 

In the dark matter potential well, the protogalactic 
gas probably settles into a two-phase medium consisting 
of (1) a hot and tenuous gas at the virial temperature, 
That — 1.7 x 10 6 K and (2) a collection of much colder self- 
gravitating clouds, pressure-bound by the hot medium in 
which they are embedded. These cold and dense clouds are 
often considered as the formation sites of the halo stars and 
GCs. In fact, this description of the protoGalaxy is encoun- 
tered in widely different pictures for halo GC formation. For 
instance, in the frame of the Galaxy formation model relying 
on a monolithic collapse, Fall & Rees (1985) suggested that 
the development of such a two-phase medium is promoted 
by a thermal instability. On the other hand, adopting the 
hierarchical picture of Galaxy formation as the framework 
of their model, Harris & Pudritz (1994) proposed that the 
formation of GCs took place in the densest parts of large 
protogalactic fragments, of mass 10 8 — 10 9 Mq (i.e., their Su- 
per Giant Molecular Clouds, SGMC), these SGMCs being, 
as in the Fall & Rees (1985) model, embedded in a hot back- 
ground at the virial temperature. 

We thus assume that the protogalaxy includes three 
components: a dark matter corona and a set of star forming 
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Figure 2. Fit of the model p co id(D) (equation 8) to the presently 
observed mass density profile of the OH subsystem. The best (in 
the x 2 sense) fit is obtained when a downwards scaling of 2 X 10~ 4 
is applied to p co i d (D) 
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Figure 3. Comparison between the radial mass density profile of 
the OH subsystem and the model mass density p co ld(D), consid- 
ering different values for the coefficients in equation (8) 



clouds embedded in a hot and diffuse background, and we 
write: 

POai(D) = p co i d (D) + p hot {D) + pdm(D) . (3) 

In this equation, pcai, pcoid, phot and p DM are the mass 
densities of the protoGalaxy, the cold phase, the hot phase 
and the dark matter, respectively. 

The density profile pGai of the protoGalaxy is conve- 
niently described as a singular isothermal sphere: 

V 2 1 

pGal = 4^GD^' (4) 
where V c is the circular velocity of the gas in the dark mat- 
ter potential well of the Galaxy and G is the gravitational 
constant. Since we are mainly interested in the proto-Milky 
Way, we adopt V c =220 km.s -1 . 

Although physically unmotivated, we adopt the usual 
description of the spatial distribution of the dark matter 
mass in the Galaxy, that is: 

PODM 



Pdm(D) 



1 + 



(5) 



As for the central density and the softening length, we adopt 
the values of the Caldwell & Ostriker (1981) model: podm = 
13.72 x lO 6 M .kpc- 3 and D DM = 7.8151 kpc. 

Regarding the hot gas confining the protoglobular 
clouds, its pressure profile is expected to scale as D~ 2 (e.g., 
Murray & Lin 1992, Harris & Pudritz 1994): 



p mi 1^5x10-9 _ 2 
Phot(D) = dyne.cm , 

kpc 



(6) 



where the coefficient comes from Murray & Lin (1992). 

Combining equation (6) with the virial temperature 
Thot, we obtain the mass density profile of the hot tenuous 
gas: 



(7) 



i-r-i-. pOhot 78.38 6,, i -3 

Phot(D) = = -^2— x 10 M Q .kpc . 

kpc 

The combination of equations (3), (4), (5) and (7) pro- 
vides us with an estimate of the radial mass density profile 



of the cold gas, that is, the gas out of which halo GCs pre- 
sumably formed: 



, , , 890.7 78.4 

Pcold(D) = | — — - 



13.7 



D 2 



xlO 6 M .kpc~ 3 .(8) 



Figure 2 displays equation (8), as well as the observed 
mass density profile of the OH GCS. The best (in a least- 
squares sense) agreement with the observed profile is ob- 
tained when the density of the protogalactic cold gas is 
scaled downwards by a factor 2 x 10~ 4 . x 2 and Q values 
for the two bin sizes are provided in Table 3 (case a). Ob- 
viously, the shape of p co id(D) provides a good fit to the 
observed distribution of the OH mass in our Galaxy. The 
factor 2 x 10~ 4 thus represents the formation efficiency of 
the halo bound stellar clusters which have managed to sur- 
vive a Hubble time in the tidal fields of the Milky Way. 
This factor of course represents a lower bound to the GC 
primordial formation efficiency. As already highlighted by 
McLaughlin (1999), the next Section confirms that this one 
is unlikely to have been exceedingly larger. 

Additionally, we note that the model, as described by 
equation (8), fits nicely the extent of the OH, as the cold 
gas density drops sharply at D ~ 40 kpc. The outer limit 
of the OH would thus arise from the depletion of cold gas 
at this galactocentric distance. However, we note that this 
agreement may be a mere coincidence only as it heavily de- 
pends on the choice of the parameters pohot, Podm and Ddm 
in equation (8). Figure 3 illustrates how variations in the 
adopted values of the parameters affect the shape of the cold 
gas profile. The dashed-dotted and dotted curves represent 
equation (8) with the hot gas coefficient pohot increased and 
reduced by 50 per cent, respectively, while the dashed curve 
shows pcoid(D) if the dark matter coefficients of Caldwell & 
Ostriker (1981) are replaced by podm — 30 x 10 6 M .kpc~ 3 
and Ddm = 5 kpc (e.g., Baumgardt 1998). In the last two 
cases, the gas density decreases smoothly with galactocen- 
tric distance, without any sharp drop in the gas density. In 
that case, the outer limit of the OH might characterize a gas 
density threshold below which the formation of bound clus- 
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Figure 4. Comparison between the same set of curves as in the 
previous figure and the spatial distribution of the OH subsystem 
modified in the following way: exclusion from the initial sample of 
the GCs with extended structures (filled points) and of the GCs 
with kinematics more typical of the YH (open points) 



ters is either inhibited or delayed. In the first case, the gas 
density drops sharply at closer galactocentric distance (i.e., 
D ~ 25 kpc). Nevertheless, we emphasize that all four curves 
fit well the observed OH profile in the range 1-20 kpc, thus 
accounting for the vast majority of the OH GCs (less than 
ten per cent of the OH GCs are located beyond D ~ 20 kpc) . 
This strengthens our hypothesis following which the mass 
density profile of the OH GCS is tracing the cold baryonic 
material available to the star formation process some 15 Gyr 
ago. 

An additional issue of potential concern is that a small, 
but non negligeable, fraction of the OH may itself be an ac- 
creted component. Several of the GCs in nearby dwarf galax- 
ies have colour-magnitude diagram characteristics that are 
indistinguishable from OH GCs (see, e.g., Grebel & Gal- 
lagher 2004 and references therein). Furthermore, Mackey 
& Gilmore (2004) emphasize that, although the overwhelm- 
ing contribution of clusters from external galaxies would be 
to the Galactic YH subsystem, there may also be a non- 
zero contribution to the OH subsystem, which they esti- 
mate of order 15 per cent or, equivalently, a dozen of GCs. 
Such accreted OH cluster candidates may be identified on 
the basis of either an extended core radius or spatial mo- 
tions more typical of YH objects, although these alone are 
not definite indicators of an extra-Galactic origin. Mackey 
& Gilmore (2004) find five OH GCs in the first category 
(NGC 6809, 6101, 7492, 5897 and Pal 15) and 6 GCs in 
the second (NGC 1904, 2298,5024,5904,6205,7089). We have 
thus refitted equation (8) to the observed OH distribution, 
after rejecting (1) the five GCs with extended structure, (2) 
all eleven GCs quoted above plus Pal 2 (see below). These 
two observed spatial distributions are illustrated in Fig. 4. 
In the first case (filled symbols) , one hardly detects a differ- 
ence with respect to the distribution of the whole OH GCS 
(as in Figs. 1-3). Accordingly, the parameters of the fit are 
much similar to previously (compare case (a) and case (b) 
in Table 3). In a second step (open symbols in Fig. 4 and 
case (c) in Table 3), we have removed the complete sample 



Table 3. Results of fitting the spatial distribution of the baryonic 
cold matter (equation 8) on the observed mass density profile of 
the OH subsystem: (a) all the GCs defined as OH objects on the 
basis of their location in a [Fc/H] vs HBR diagram (Table 1 and 
Figs. 1, 2 and 3); (b) same as (a) but excluding five GCs with 
extended structures (filled symbols in Fig. 4); (c) same as (a) but 
excluding objects with extended structure or extreme kinematics 
(open symbols in Fig. 4) (see text for details) 





ALogio-D 


ALog 1() p 


x 2 


Q 


(a) 


0.1 


-3.69 ± 0.06 


15.0 


0.45 




0.2 


-3.70 ± 0.08 


6.3 


0.51 


(b) 


0.1 


-3.70 ± 0.06 


14.5 


0.49 




0.2 


-3.71 ± 0.08 


6.3 


0.50 


(c) 


0.1 


-3.78 ± 0.07 


20.6 


0.11 




0.2 


-3.81 ± 0.08 


11.1 


0.14 



of 11 OH GCs suspected of having been accreted as well as 
Pal 2. In the [Fe/H] vs HBR diagram, this one is located 
at the frontier between the OH and the YH groups and its 
nature is thus ill-determined. While we have previously con- 
sidered it as an OH GC, we now sort it in the YH group. 
Removing the so-defined sample of 12 GCs, the changes are 
more significant but still, the modified observed distribution 
is satisfactorily fitted by equation (8), the main point being 
that the model distribution overestimates the mass density 
of the OH subsystem in the region D ~ 20 kpc. We thus 
conclude that our fit is robust, even though the actual ra- 
dial distribution of the OH mass remains slightly uncertain 
as a few old clusters may be accreted objects and, thus, in- 
terlopers with respect to the genuine initial Galactic GCS. 

In what follows, we consider the OH mass distribution 
displayed in Figs. 1 to 3 and equation (8) as our fiducial ob- 
served and theoretical radial mass density profiles, respec- 
tively. 



4 EVOLUTION OF THE SPATIAL 

DISTRIBUTION OF THE OH MASS 

As announced in Section 2, we now put on a firmer foot 
the hypothesis following which the initial GCS mass density 
profile is reasonably approximated by what it presently is. 

Evolutionary processes act more efficiently upon low- 
mass GCs as well as on GCs located in the inner Galactic 
regions. In other words, the evolution with time of the GCS 
has been mostly determined by the initial spatial distribu- 
tion of the clusters in the Galactic halo as well as by their 
initial mass spectrum. In this respect, it is interesting to 
note that the sharp contrast between the fraction of surviv- 
ing clusters and the ratio of the final to the initial total mass 
in clusters as obtained by McLaughlin (1999) (see Section 
2) arises from the choice of a power-law mass spectrum with 
a steep (i.e., -2) slope and probing down to very low-mass 
(i.e., of order 10 3 M ). Other mass spectra could thus lead 
to different results. 

How the initial mass spectrum of the Galactic halo GCs 
looked like is a much debated issue. As mentioned earlier, the 
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Table 4. Results of our simulations for different GC initial mass functions (IMF, see text for details). Fj\r and Fm are the fraction of 
surviving GCs and the ratio of the final to the initial mass in GCs, respectively. Also given are the results of fitting the radial distributions 
evolved during 15 Gyr to the data for the two different binnings considered, that is, the vertical logarithmic shift ALog 10 p required to 
obtain the best match between the model and the data and the corresponding \ 2 and Q(z//2, x 2 /2) values 



Run 


GC IMF 


Fat 


F M 




ALog 10 D=0.1 






ALog 10 D=0.2 












ALog 10 p 


x 2 


Q 


ALog 10 p 


x 2 


Q 


[1] 


Gaussian (Mean=5.03; <r=0.66) 


0.75 


0.66 


-2.11 


20.7 


0.15 


-2.14 


7.1 


0.40 


[2] 


Power-law (m iotl> =10 3 M©) 


0.06 


0.39 


-0.54 


27.3 


0.03 


-0.56 


10.9 


0.14 


[3] 


Power-law (m; OTO =10 4 M© ) 


0.46 


0.51 


-1.44 


27.2 


0.03 


-1.46 


10.9 


0.14 


[4] 


Power-law (m; otl) =10 5 M©) 


0.96 


0.68 


-2.30 


20.4 


0.16 


-2.31 


6.9 


0.44 



observed luminosity spectrum dN/dL of numerous cluster 
systems recently formed in starburst and merging galaxies 
is a power-law. However, whether the luminosity spectrum 
dN/dL constitutes a faithfull mirror of the underlying mass 
spectrum dN/dm is questionable. Both shapes agree only 
if the variations of the mass-to-light ratio from cluster to 
cluster are limited (as is roughly the case for the Galactic 
halo GCS). Such a requirement may not be true for ongoing 
or recent starbursts: the formation duration of such systems 
may be a significant fraction of the system's median age and, 
thus, age spread effects among the young star cluster popu- 
lation may not be negligible. Being an age related quantity, 
the mass-to- light ratio can no longer be considered as a con- 
stant and the shape of the luminosity spectrum may differ 
substantially from the shape of the mass spectrum. Raising 
this issue, Meurer et al. (1995) and Fritze v. Alvensleben 
(1998, 1999) showed that systems of young GCs could dis- 
play a power-law luminosity spectrum while the underly- 
ing mass spectrum is a broken power-law or, equivalently, a 
gaussian when the binning is logarithmic (dAVdlogm). The 
reason for this is that, owing to the fading with time of the 
cluster luminosity, high mass clusters can be observed over 
a wide range of ages, while low-mass ones are detectable at 
young ages only. As a result, low-mass clusters are under- 
represented in the observed mass spectrum. Thus, this one 
can be described by a two-index power-law (or equivalently 
a bell-shaped mass function) as it does not raise in the low- 
mass regime as steeply as in the high-mass regime. 

Vesperini (1998) has indeed shown that the halo GCS 
could have started with a gaussian initial mass function. 
Building on N-body simulations performed by Vesperini & 
Heggie (1997), he demonstrated the existence of a quasi- 
cquilibrium GC mass function, that is, a mass function 
whose initial gaussian shape and parameters (mean and 
standard deviation) are preserved during the entire evolu- 
tion through a subtle balance between disruption of clus- 
ters and evolution of the masses of those which survived, 
even though a significant fraction of the GCs is destroyed. 
Interestingly, the mean and the standard deviation of this 
gaussian equilibrium mass function (log (m/M©) = 5.03 and 
a — 0.66) are remarkably close to those of the mass function 
of the halo GCS. Thus, the initial GC mass function may 
have been very similar to what it is today. Obviously, the 
study of the temporal evolution of the GC mass spectrum 
is not enough to unveil its initial shape as both a power-law 



and a gaussian mass functions evolve into the current GC 
mass function after a Hubble time. 

If we assume that the halo GC initial mass spectrum 
was similar to what is observed today in starbursts and 
mergers, to derive the latter, one should estimate the in- 
dividual ages for all star clusters. This can be done by, for 
instance, comparing the observed broad-band photometry to 
the colours generated by spectral evolution synthesis mod- 
els (see, e.g., Parmentier, de Grijs & Gilmore 2003). These 
models will then provide an estimate of the mass-to-light 
ratio for each cluster as a function of age, eventually en- 
abling one to derive the intrinsic mass spectrum underlying 
the observed luminosity spectrum. Unfortunately, even such 
studies lead to unconclusive results. Analysing the system 
of young GCs in the Antenna merger NGC 4038/39, Zhang 
& Fall (2001) show that the very young clusters (i.e., with 
age less than 150 Myr), those which have remained unaf- 
fected by the various dynamical effects, are distributed in 
mass according to a pure power-law with slope —2. On the 
other hand, de Grijs et al. (2003) derive a gaussian mass 
function dAVdlogm (or equivalently a two-index power-law 
when binning the data in m) for an equally young system 
in the nearby starburst galaxy NGC 3310. 

Therefore, at the present stage, neither the theory of 
the dynamical evolution of a GCS, nor the observations of 
young GCs in starburst galaxies can help distinguishing be- 
tween a power-law or a gaussian mass function for the GC 
initial population. Considering the very unclear issue of this 
debate, we have computed the temporal evolution of the 
spatial distribution of the GCS mass around the Galactic 
centre for various mass functions. 

To evolve the spatial distribution of the mass of a GCS 
from the time of its formation up to an age of 15 Gyr, we 
adopt the analytic formulae of Vesperini & Heggie (1997) 
which supply at any time t the mass m of a GC with initial 
mass rtii and moving along a circular orbit at a galactocen- 
tric distance D. These relations have been obtained by fit- 
ting the results of a large set of N-body simulations in which 
Vesperini & Heggie (1997) take into account the effects of 
stellar evolution as well as two-body relaxation, which leads 
to evaporation through the GC tidal boundary. Disc shock- 
ing can also be included (see below). In order to take into 
account dynamical friction, GCs whose time-scale of orbital 
decay (see, e.g., Binney & Tremaine 1987) is smaller than t 
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Figure 5. Evolution with time of the radial mass density profile 
of a GCS whose mass is initially spatially distributed according to 
equation (8) (dashed curve). Results with (dashed-dotted curve) 
and without (plain curve) disc-shocking (DS) are shown for three 
different initial GC mass spectra. The plain curve is vertically 
shifted to match the observed data (see the quantity ALog 10 p in 
Table 4) 



are removed from the GCS at that time. A summary of the 
method is provided in Vesperini (1998, his Section 2). 

The temporal evolution of the mass of a GC, modelled 
as a multi-mass King model with an initial dimensionless 
concentration parameter Wo = 7 and orbiting at constant 
galactocentric distance D, is assumed to follow: 



time t is expressed in units of 1 Myr and F cm , a quantity 
proportional to the initial relaxation time, is defined as: 



p 

- - ... : 1 : p 

- 


oh(D) 
oh(D) 

1 
1 


A(LogD) = 0.: 
A(LogD) = 0.5 

WD) : 

) end (D)-noDS 
'end(D)- DS 


■ 






























































■ 




|Run 4| Power law MSp: 


, = -1.9 




,= 1E5 


M \ 



m(t) 



= 1 - 



Arris 



0.828 



(9) 



mi x D 
kiN 



(10) 



Am si! „/mi is the fraction of cluster mass lost due to 
stellar evolution (18 per cent in this particular model). The 



where rm and D are in units of 1 Mq and 1 kpc, re- 
spectively, and N is the initial number of stars in the GC. 
To take into account disc shocking, the factor O.S2S/F cw is 
merely replaced by A as defined by equation (3) of Vesperini 
(1998). 

We have distributed 20,000 GCs following a radial dis- 
tribution obeying our equation (8) and various mass spectra: 

a) a gaussian mass function diV/dlogm with parameters 
equal to those of the equilibrium mass function of Vesperini 
(1998), 

b) three power- law mass spectra AN /Am with a single 
slope of -1.9 and different lower mass limits, namely, 1E3, 
1E4 and 1E5M . 

As for the last cut-off value, Fall & Zhang (2001) indeed 
showed that the GC mass spectrum might have started with 
a truncation at mass of order 1E5 Mq , the low-mass tail of 
the currently observed GC mass distribution being formed 
as a result of the evaporation of the massive GCs located at 
short distance from the Galactic centre. 

Each panel of Fig. 5 displays the results of one of these 
simulations (the case of a power-law truncated at 1E4 M© is 
not represented as it is highly similar to the 1E3M Q case). 
The dashed, plain and dashed-dotted lines are, respectively, 
the initial mass density profile as described by equation (8), 
the radial mass density profile after a Hubble time long evo- 
lution without disc shocking (obtained through equation 9) 
and with disc shocking. The final profile without disc shock- 
ing (plain curve) has been vertically shifted to provide the 
best agreement, in the least-squares sense, with the observed 
spatial distribution of the OH mass. The amplitude ALog 10 p 
of the shift, the \ 2 an d Q values of each fit, as well as the 
fraction of surviving clusters (Fjv) and the ratio of the final 
to the initial total mass in GCs (Fm), are listed in Table 4. 
Fjv exhibits an extremely large range of variation, with more 
than 90 per cent of the clusters being disrupted if the initial 
mass spectrum is a single power-law going down to 1000 M , 
while a cut-off at 10 5 M leads to the preservation of most of 
the GCs. This large scatter in Fjv contrasts markedly with 
the limited range in Fm, as this one varies by less than a fac- 
tor 2. Therefore, considering our ignorance of the shape of 
the GC initial distribution in mass, mass-related quantities 
appear as indicators of the initial conditions more reliable 
than number-related quantities. In the frame of the initial 
density profile described by equation (8), we also note that 
the initial mass of the OH GCS was at most 2.5 times larger 
than what it currently is, thus suggesting that the bound 
cluster formation efficiency in the Galactic halo was at most 
of order 2.5 x (2 x 10~ 4 ). This result also agrees with previous 
studies (e.g., Baumgardt 1998, Vesperini 1998, McLaughlin 
1999) following which the disruption of GCs cannot be a ma- 
jor contributor to the build-up of the stellar halo. Indeed, 
assuming a mass-to-light ratio of 2.35 (see Section 2), the 
mass of the OH GC subsystem is ~ 2 x 10 7 Mq, about two 
orders of magnitude less massive than the stellar halo (~ 
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10 9 M Q , Freeman & Bland-Hawthorn 2002). 

As already quoted, the shapes of the initial and final 
mass profiles are very similar, at least beyond a galacto- 
centric distance of ~ 2kpc. It may be worth remembering 
that about one fifth of the cluster mass decrease is due to 
mass-loss associated with stellar evolution, which is inde- 
pendent of galactocentric distance. As a result, the shape 
of the radial profile is left unaffected by this gaseous mass 
loss taking place during the very early evolution of the GCs. 
The initial shape of the model profile being only moderatly 
affected by dynamical evolution, it is not surprising that the 
observed OH profile (indicated by the filled squares or cir- 
cles in Fig. 5, depending on the size of the bins) and the 
outcome of our simulations (plain curve in Fig. 5) are in 
satisfactory agreement as indicated by the Q values listed in 
Table 4. The agreement is weaker than in Table 3 however, 
especially for the power-laws with a truncation at 1E3 Mq 
and 1E4M©, but this effect is mostly driven by one bin 
only, at logD = 0.25. Should it be ignored, the agreement 
between the model and the observations would jump up to 
Q ~ 0.3 in the least favourable cases. This kick in the ob- 
served density profile arises from the presence, at almost the 
same galactocentric distance, of two of the handful of very 
massive (M v ~ -9.2) GCs, i.e., NGC 6266 and NGC 6273 
at D=1.7 and 1.6 kpc, respectively. 

In case of disc shocking (dashed-dotted curves in Fig. 5), 
the ratios Fjv and Fm are almost unaffected, being reduced 
by 4 per cent at most with respect to the simulations with- 
out disc shocking (see Table 4). The shape of the density 
profile is also preserved except at galactocentric distances 
shorter than 2 kpc where the predicted cluster density is sig- 
nificantly affected compared to the case without disc shock- 
ing. In these regions, the disc shocking model underestimates 
the observed mass density. This discrepancy may be solved if 
these very inner halo regions are contaminated by the metal- 
poor tail of the bulge GCS. Alternatively, it may be that the 
high-density high-pressure environment of these very inner 
Galactic regions promoted the formation of GCs with an ef- 
ficiency locally higher than in the overall halo. In such case, 
the original GCS mass density profile around the Galactic 
centre may have been steeper than the density profile p co u- 

It is interesting to compare the values Fn and Fm listed 
in Table 4 with those obtained in previous studies. Consid- 
ering the case of the equilibrium GC mass function (run 
[1] in Table 4), Vesperini (1998) derived Fjv ~ 0.55 and 
Fm — 0.40. As for a power- law mass spectrum with a slope 
-2 and a low-mass cut-off ~1000M© (run [2] in Table 3), 
Baumgardt (1998) indicates that the fraction of surviving 
clusters may be of order 1 per cent, while the final mass in 
GCs would represent 15 to 30 per cent of the initial mass 
content. These results, although in rough agreement with 
ours, are systematically lower. This is not surprising how- 
ever as these studies assume an initial radial distribution 
scaling as D~ 3,5 (even D -4 ' 5 for some of Baumgardt's runs), 
thus significantly steeper in the inner regions than our radial 
profile p co id (see Fig. 2) . Therefore, with respect to our sim- 
ulations, a larger fraction of the initial population of GCs is 
put on orbits closer to the Galactocentric centre where they 
are more efficiently evaporated and eventually disrupted. 



5 DISCUSSION AND SUMMARY 

We have proposed a model accounting for the spatial distri- 
bution of the mass of the OH GCS. We suggest that this GC 
subsystem formed out of the cold baryonic gas present in the 
Galaxy some 15 Gyr ago, that is the gas leftover once the 
contribution of the dark matter is removed, with a roughly 
constant formation efficiency throughout the halo. In or- 
der to test our hypothesis, we have built the presently ob- 
served radial distribution of the OH GCS, making use of 
the latest classification of the halo GCs between the OH 
and YH subsystems (Harris 1996, updated 2003, Mackey & 
Gilmore 2004). We have found good agreement between our 
suggested model (equation 8) and the observed OH mass 
distribution (see Fig. 2 and Table 3). The fit is robust with 
respect to uncertainties in the observed distribution of the 
mass (i.e., depending on whether a few OH GCs may be ac- 
creted objects as the YH GCs) as well as with respect to the 
uncertainties in the parameters involved in equation (8) . 

We have made one more step, demonstrating that the 
GCS mass density profile, as described by equation (8), has 
not been severly affected by a Hubble time-long evolution 
in the tidal fields of the Milky Way, at least beyond galac- 
tocentric distances of order 2 to 3 kpc. This result is only 
weakly dependent on the shape of the initial GC mass spec- 
trum, which is still very poorly known. In fact, numerous 
studies found that both an initial gaussian mass function 
(i.e., a two-index power law mass spectrum) and an initial 
power-law mass spectrum will evolve into the presently ob- 
served gaussian mass function of GCs. Therefore, to recover 
the faint end of the original mass function of GCs on the 
basis of its temporal evolution only is not feasible. We have 
thus considered various mass functions in order to demon- 
strate the robustness of the GCS mass density profile with 
respect to the dynamical evolution. 

Our simulations make use of the analytic formulae ob- 
tained by Vesperini & Heggie (1997), and the starting point 
consists in assuming a given radial distribution and a given 
mass function for the GCs. Doing so, we neglect finer de- 
tails such as the non-circularity of the GC orbits as well as 
the initial concentration of GCs. As for the first point how- 
ever, we recall that our main interest in this paper is the 
OH subsystem, which shows less extreme kinematics than 
the YH GCs (see, e.g., Dinescu, Girard & van Altena 1998). 
Thus, while the assumption of circular orbits for the GCs 
is clearly a simplifying one, this issue is likely to be less 
critical than if we have dealt with the whole halo GCS or 
with the YH only. Regarding the second point, we note that 
the equation (9), which describes the temporal decrease of 
the cluster mass according to Vesperini & Heggie (1997), 
implicitly assumes an initial concentration Wo = 7 for all 
GCs. The existence of GCs with lower initial concentration 
is thus clearly neglected. Such low concentration GCs may 
have been initially present and destroyed relatively quickly 
by a combination of stellar winds and tidal limitation with 
little or no help from two-body relaxation or gravitational 
shocks. In the inner Galactic regions (say within the Solar 
Circle), that is, where most of the GCS mass is located, even 
high-mass GCs are unable to survive if their initial concen- 
tration is low (Vesperini 1997). By assuming that all GCs 
start with a concentration Wo = 7, we may thus overesti- 
mate the survival capacity of some high mass GCs and, thus, 
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underestimate the decrease with time of poh(D). However, 
it is worth pointing out that GCs exhibit a correlation be- 
tween their luminous mass estimate and their concentration 
in the sense that more massive GCs are more concentrated 
(van den Bergh 1994). This correlation is likely to be of 
primordial origin as it is stronger at larger galactocentric 
distances (Bellazini et al. 1996), that is, where evolutionary 
processes act on longer time-scales and, thus, where mem- 
ory of the initial conditions is better preserved. Moreover, 
simulations show that high-mass low-concentration GCs are 
destroyed in the inner halo but manage to survive at large 
galactocentric distances (i.e., D > lOkpc, Vesperini 1997). 
As they are not observed even beyond the Solar Circle, high- 
mass low-concentration GCs were probably not part of the 
initial GCS and, therefore, to neglect them should not affect 
the results displayed in Fig. 5 and Table 4. 

Bearing these caveats in mind, we have shown that the 
shape of the evolved radial mass density profile still nicely 
mirrors the initial one, even though a significant fraction 
of the GCs has been destroyed during a Hubble time long 
evolution. Not surprisingly thus, the observed spatial distri- 
bution of the OH mass and the distributions arising out of 
our simulations are also found to be in good agreement (see 
Fig. 5 and Table 4). This strengthens our hypothesis fol- 
lowing which the distribution of the gas GCs formed from 
obeys equation (8) and, also, that the dark matter, what- 
ever its nature, was not involved in the earliest stages of star 
formation. 

This good agreement also suggests that the flattening 
of the observed distribution of GCs around the Galactic cen- 
tre may not result from tidal disruption processes only but 
may be, at least partly, of primordial origin. Indeed, in the 
region 1 to 3kpc, our suggested cold matter profile shows 
a flattening (slope ~ —2) which contrasts with the over- 
all shape, reasonably approximated by a slope of ~ —3.5 
(see Fig. 2). The situation is still unclear however as, at an 
age of 15 Gyr, our mass distribution underestimates the ob- 
served mass density profile of the GCS for D < 2kpc, espe- 
cially when disc shocking is included in the simulations (see 
Fig. 5). Possible solutions of this discrepancy include the 
contamination of the observed distribution by the metal- 
poor tail of the bulge/disc system and/or star formation 
efficiencies higher in the inner regions of the Galaxy owing 
to higher pressures and densities. As a result of this latter 
suggestion, the initial mass distribution of the GCs in the 
Galactic central regions may have been steeper than indi- 
cated by equation (8). Apart from these uncertainties, the 
satisfactory agreement between the model and the observed 
radial density profile of the OH subsystem, beyond D ~ 2 
to 3kpc, suggests that the formation of GCs took place with 
a formation efficiency almost constant throughout the OH. 

This paper has been typeset from a TfrjX/ F/TfrjX file prepared 
by the author. 
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